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ABSTRACT 

We present Keck/DEIMOS spectroscopy and CFHT/MegaCam photometry for the Milky Way 
globular cluster Palomar 13. We triple the number of spectroscopically confirmed members, including 
many repeat velocity measurements. Palomar 13 is the only known globular cluster with possible 
evidence for dark matter, based on a Keck/HIRES 21 star velocity dispersion of a = 2.2 ±0.4 km s~^. 
We reproduce this measurement, but demonstrate that it is inflated by unresolved binary stars. For 
our sample of 61 stars, the velocity dispersion is cr = 0.7j^Q gkm s^^. Combining our DEIMOS data 
with literature values, our final velocity dispersion is cr = 0.4^q'3 km s^^. We determine a spectroscopic 
metallicity of [Fe/H]= —1.6 ± 0.1, placing a 1-sigma upper limit of cr[Fc/H] ^ 0.2 dex on any internal 
metallicity spread. We determine Palomar 13's total luminosity to be My = —2.8 ± 0.4, making it 
among the least luminous known globular clusters. The photometric isophotes are regular out to the 
half-light radius and mildly irregular outside this radius. The outer surface brightness profile slope 
is shallower than typical globular clusters (S oc r'',?7 = — 2.8±0.3). Thus at large radius, tidal 
debris is likely affecting the appearance of Palomar 13. Combining our luminosity with the intrinsic 
velocity dispersion, we find a dynamical mass of of A/1/2 = l-Sii'a ^ lO^Af0 and a mass-to-light ratio 

of M/Lv — 2.4^2 4-^0/^0- Within our measurement errors, the mass-to-light ratio agrees with the 
theoretical predictions for a single stellar population. We conclude that, while there is some evidence 
for tidal stripping at large radius, the dynamical mass of Palomar 13 is consistent with its stellar 
mass and neither significant dark matter, nor extreme tidal heating, is required to explain the cluster 
dynamics. 

Subject headings: globular clusters: individual (Palomar 13) - galaxies: kinematics and dynamics 
-dark matter 



1. INTRODUCTION 

Globular clusters are among the oldest structures in 
the universe. Current models favor globular cluster for- 
mation via the collapse of cold gas withi n a larger galac- 
tic environment (1^ 11 & Recs 1985; Kra vtsov fc GnedinI 
l2005HMuratov fc Gnedini201Qil . Recently, improved ev- 
idence for internal abundance ratios of elements heav- 
ier than hydrogen i n some of the most luminous glob- 
ular clusters (e.g iLee et all 120091: iCohen et all 120101: 
iMilone et al.ll201ll ) has renewed interest in the idea that 
these objects might instead have formed at the center 
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of their own dark matter halos (iPeeblesI I1984D or are 
the s tripped nuclei of dwarf galaxies ( Bekki &: FreemanI 
|2003| ) , thus providing a deeper gravitational potential for 
the retention of enriched material. Observational evi- 
dence, however, argues against dark matter in globular 
clusters, at least at t he present time. Ph otometric stud- 
ies of cluster profiles (Conroy et al."2010) and dynamical 
studi es of Milky Way g lobular clusters (Baun igardt et al\ 
[20091: iLane et al.ll2010[lHankev fc Colell201ir show no ev- 
idence for dark matter, at least out to the tidal radius in 
these systems. 

An interesting possible exception to the conclusion 
that globular clusters do not contain dark matter is Palo- 
mar 13 (Pal 13), a lo w luminosity Milky Way globular 

cluster (My 3.7: iHarrid IMol ) at a Galactocentric 

distance of 25 kpc. The combined size and luminosity of 
Pal 13 (see Figure [Ij is more similar to the Milky Way 
ultra-fai nt galaxies than a typical Milky Way globular 
cluster (|Martin et a l. 2008). Motivat ed in part by its 
low luminosity and relatively large size, ICote et al.l(l200l 
hereafter C02) measured velocities for 21 member stars 
in Palis based on Keck/HIRES spectroscopy. Their ve- 
locity dispersion of 2.2 ± 0.4 km s~^ is well in excess of 
that predicted from the stellar mass of Pal 13 alone. In 
contrast, velocity dispersions of six other extended outer 
halo globular clusters measured by these authors are con - 
sistent with their stellar masses (iBaumgardt et aLll2009f) . 
C02 also measure a shallow surface brightness profile for 
Pal 13 and a King tidal radius that is significantly larger 
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than the predicted tidal radius based on the stellar mass 
of the cluster. C02 interpreted the inflated velocity dis- 
persion and shallow luminosity density profile of Pal 13 as 
due either to tidal heating during a recent perigalacticon 
passage, implying that the cluster is not in dynamical 
equilibrium, or the presence of dark matter. If Pal 13 is 
in dynamical equilibrium, its measured velocity disper- 
sion implies the highest mass-to-light ratio known in a 
globular cluster, M/Ly = Mq/ Lq. 

The proper motion of Pal 13, measured bv' Siegel et al.l 
(pbOli) via photographic plates, implies that this cluster 
is on a highly eccentric orbit near perigalacticon. This 
is consistent with the interpretation that Pal 13 is in the 
final stages of tidal disruption, and its increased veloc- 
ity dispersion is due to tidal heating. However, N-body 
simulations by Kiipper ct al. (2011) are unable to repro- 
duce either the inflated velocity dispersion or shallow 
density profile using the measured orbit. Interestingly, 
these authors can reproduce the photometric properties 
of Pal 13 assuming that instead it is near apogalacticon, 
i.e., at the farthest point in its orbit. In this case, the 
observed surface brightness profile is not due to tidal 
heating, but rather due to projection effects as unbound 
stars, removed in the course of slow tidal evaporation, are 
compressed near apogalacticon. In these simulations, bi- 
naries are still required to explain the observed velocity 
dispersion. 

A high unresolved binary star fraction is a third pro- 
posed explanati on for the large ve locity dispersion ob- 
served in Pal 13. iClark et al.l (|2004[ ) estimated the Pal 13 
binary star fraction as 30 ± 4% based on a secondary 
red main sequence and the presence of a large number 
of blue straggler stars. This is higher than observed for 
more massive globular clusters, but consis tent with that 
of lower mass globular clusters. iBlecha e~al. (2004 ) mea- 
sured a lower velocity dispersion (0.6 — 0.9 ± 0.3 km s~^) 
for Pal 13 based on eight member stars. Seven of these 
stars overlap with C02 (note Blecha et al. state only six), 
including three stars which are identified as velocity vari- 
ables in one or both samples. Given the small number of 
stars, it is yet unclear whether binaries can indeed pro- 
vide an explanation for the high velocity dispersion in 
Pal 13 without invoking tidal heating or the presence of 
dark matter. 

Here we present Keck/DEIMOS spectroscopy and 
CFHT/Megacam imaging of the Milky Way globular 
cluster Pal 13. The paper is organized as follows: in §[2] 
we discuss data acquisition and reduction for both our 
spectroscopic and photometric datasets, including esti- 
mates of structural parameters. In §[3] we discuss the 
techniques used to identify spectroscopic members. In 
§13] we determine the velocity dispersion, mass and mass- 
to-light ratio of Pal 13, and compare to previous results. 
In §|4] we also examine the blue straggler population of 
Pal 13. Finally, in §[5] we interpret our results, conclud- 
ing that the kinematics of Pal 13 are consistent with its 
stellar mass and find only mild evidence for tidal effects 
at large radii. 

Throughout this paper, we adopt a distance to Pal 13 
of 24.3 kpc (Rgc = 25.3 kp c: iCote et al. 200^ and red- 
dening of E(B-V) = 0.11 (jSchlegel et al...l998h . These 
values are in good agreement with our own isochrone fit- 
ting of the photometry described in § 12.2.11 



2. DATA ANALYSIS 
2.1. Spectroscopic Targeting 

We select spectroscopic targets in Pal 13 based on 
Canada-France-Hawan Telescope (CFHT) CFH12K and 
Keck/LRIS imaging as described in C02. We supple- 
ment these data at large radius with wider field pho- 
tometry from P. Stetson's online databascFl. Candi- 
date Pal 13 member stars were identified in V- and I- 
bands, and V- and i?-band when /-band was not avail- 
able. Target s were selected based on proximity to a 
iGirardi et al.l ([2002) theoretical isochrone for an age of 
12Gyr and a metalhc i ty of [Fe/H]= - 1.7 (jZinn fc Diad 
[T98a iFriel eTallfTQSl iCote et al.ll200l . This isochrone 
was used for target selection only; a slightly better fitting 
isochrone is used in the subsequent analysis. We deter- 
mine a star's distance to the isochrone via the quantity 
d = ^[{V - /), - (V^ - /)iso]2 + (K - V^^oY. stars on 
the red giant branch {V < 20.5) within d = 0.1 mag 
of the isochrone were given the highest priority for spec- 
troscopy; stars within 0.2 mag of the isochrone were given 
lower priority. On the subgiant branch and main se- 
quence {V > 20.5) and on the horizontal branch the se- 
lection region was widened to within 0.2 and 0.3 mag 
of the isochrone for the higher and lower priority cate- 
gories, respectively. Candidate blue straggler stars were 
chosen in a box blueward of and brighter than the main 
sequence turnoff /subgiant branch region. Following the 
photometric selection, spectroscopic targeting priorities 
were set to favor k nown radial velocity members from 
C02 and IBlecha et al. (2004 ), an d likely proper motion 
members from lSierel et al., (|2001l ). 

2.2. Spectroscopic Observations and Data Reduction 

The spectroscopic data were taken with th e Keck II 10- 
m tel escope and the DEIMOS spectrograph ()Faber et al.l 
|2003() . Eight multislit masks were observed in Pal 13 in 
September 2008, with a ninth mask obtained in October 
2009. Exposure times, observation dates, and other ob- 
serving details are given in Table 1. The masks were ob- 
served with the 1200 line mm~^ grating covering a wave- 
length region 6400 — 9100A. The spectral dispersion of 
this setup was 0.33A pixel^^, equivalent to R=6000 for 
our 0".7 wide slitlets, or a FWHM of 1.37A. The spa- 
tial scale was 0''12 per pixel. The minimum slit length 
was 5" which allows adequate sky subtraction; the mini- 
mum spatial separation between slit ends was 0''4 (three 
pixels) . 

Spectra were reduced using a modified version of the 
spec2d software pipeline (version 1.1.4) developed by the 
DEEP2 team at the University of California-Berkeley for 
that survey. A detailed des cription of the reductions 
can be found in iSimon fc G eha (2007). The final one- 
dimensional spectra were rebinned into logarithmic wave- 
length bins with 15 km s~^ per pixel. Radial velocities 
were measured by cross-correlating the observed science 
spectra with a series of high signal-to-noise stellar tem- 
plates. We employ the telluric absorption lines in the 
stellar spectra to calculate and apply a correction for the 
radial-velocity errors that arise from a mis-centering of 
the star in the spectrograph slit. Each science spectrum 

^'^ See |http://cadcwww.dao.nrc.ca/community/STETSON/standards7] 
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Fig. 1. — Absolute magnitude versus half-light radius for all Milky Way globular clusters (blue triangles, data from Harris 2010), dSphs 
(green circles) and ultra- faint galaxies (red circles). The open triangle represents Pal 13 based on literature values, the blue star is the 
position of Pal 13 based on quantities we determine based on deeper CFHT/MegaCam photometry. 



is cross-correlated with a hot stellar template -where the 
telluric absorption lines are the dominant spectral fea- 
tures. Since these lines are illuminated by the stellar 
flux, this defines the star's photocenter -within the slit in 
velocity space. We applied both a telluric and heliocen- 
tric correction to all velocities presented in this paper. 

We determined the random component of our velocity 
errors using a Monte Carlo bootstrap method. Noise -was 
added to each pixel in the one-dimensional science spec- 
trum, we then recalculated the velocity and telluric cor- 
rection for 1000 noise realizations. Error bars are defined 
as the square root of the variance in the recovered mean 
velocity in the Monte Carlo simulations. The system- 
atic contribution to th e velocity error "was determined by 
iSimon fc Gehal (|2007D to be 2.2 km based on repeated 
independent measurements of individual stars and we 
refer the reader to this paper for further details. The 
systematic error contribution is expected to be constant 
as the spectrograph setup and velocity cross-correlation 
routines are identical. We added the random and sys- 
tematic errors in quadrature to arrive at the final veloc- 
ity error for each science measurement. Radial velocities 
-were successfully measured for 482 of the 566 extracted 
spectra across the nine observed DEIMOS masks. The 
fitted velocities "were visually inspected to ensure relia- 
bility. We identified 90 sources as background galaxies or 
QSOs, and obtained repeat measurements for 62 stellar 
sources. Our final stellar sample consists of 392 measure- 
ments of 306 unique stars. 

2.2.1. Deep CFHT Megacam Imaging of Palomar 13 

After our spectroscopic observations, we obtained sig- 
nificantly deeper images of Pal 13 using the CFHT Mega- 
Cam imager. These images "were taken as part of a 
larger program aimed at obtaining deep "wide-field imag- 
ing of all bound stellar over-densities in the Galactic halo 
beyond 25kpc (R. Muiioz et al. 2011, in preparation). 
Pal 13 was observed on the night of UT July 22, 2009 un- 
der dark conditions with typical seeing of 0.7 — 0.9". We 
obtained six dithered exposures of 360 seconds in both 
the g- and r-bands centered on Pal 13. A standard dither- 
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Fig. 2.— CMD for CFHT photometric data within 5' of PallS's 
center. Several Girardi et al. (2002) isochrones are plotted with 
[Fe/H]= — 1.6dex for various ages shown in the legend. 



ing pattern was used to cover both the small and large 
gaps present between chips. 

Data from MegaCam were pre-processed 
prior to release using the "ELIXIR" package 
i|Magnier fc Cuillandrd [20Q1 . This pre-process m- 
cludes bad pixel masking, bias subtraction, flat fielding 
as well as preliminary photometric and astrometric 
solutions. We carried out subseq uent data re d uction 
following the method detailed in IMuhoz et al.l (|2010[ ) 
which reli es on the DA OPHOT / Allstar / ALLFRAME 
packages (jStetsonl Il994f ). To refine the astrometric 
solutions included in the iniage h eaders, we have used 
the SCAMP package (|Bertinll2006D using as reference the 
USNO-Bl catalog, which yielded astrometric residuals 
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of rms ~ 0.4". 

Photometric calibration is achieved by comparing di- 
rectly to Pal 13 data fr om the SDSS-Data Release 8 
(DR8: lA"ihara et ani2011h . To determine zero points and 
color terms we used stars with magnitudes in the range 
18 < r < 21 and -0.25 < g - r < 1.0, providing 2950 
matches between our catalog and SDSS. Final solutions 
for the zero points and color terms yield uncertainties 
of ~ 0.004 mag and ^ 0.007 mag, respectively, for both 
passbands. Our photometric catalog reaches 90% com- 
pleteness at r ^ 24.1. 

2.2.2. Improved Structural Parameters for Palomar 13 

Using our CFHT/Megacam catalog, we recalculated 
structural parameters for Pal 1 3 using a Max i mum Like- 
lihood method as described in IMuhoz et al.l ()2010D . In 
this method, a fixed analytic profile is assumed and its 
free parameters are fit using all the available stars in 
the CMD selection window. This avoids the need to bin 
or smooth data, providing more r eliable estimates fo r 
systems with low numbers of stars ()Martin et al.l [20081 ) . 
We simultaneously estimated the scale length, the cen- 
tral coordinates, the ellipticity, the position angle and 
the background de nsity of Pal 13. We assumed an em- 
pirical King profile ()KinEjll962l ). for which we calculated 
a King core (re) and tidal (rt) radius. We additionally 
fit a NFW (iNayarro et al.lll997[ ) and a Plummer profile 
ummeiill911[ ). a good representation for dwarf galax- 
ies, but found that these were poor fits to the light distri- 
bution shape of Pal 13. The best-fitting King parameters 
are = 0.42 ± 0.06' (3.0 ± 0.4 pc) and rt = 13.9 ± 1.5' 
(98 ± 11 pc). The two-dimensional half-light radius of 
this King profile is ri/2 = 1.27 ±0.16' (9.0±l.lpc). The 
overall ellipticity of Pal 13 is consistent with being zero 
and therefore its position angle is poorly constrained. 

In the left panel of Figure[3l we show the projected den- 
sity profile for Pal 13, overplotting our best-fit King pro- 
file. Error bars on the data points include uncertainty in 
determining the background level which dominates out- 
side r ~ 7'. We emphasize that the King profile shown 
is not a direct fit to the data points, but constructed 
with the best fit parameters obtained through the maxi- 
mum likelihood estimator. We compared our best fitting 
King profile to that determined in C02. If we parame- 
terize the outer slope of the surface brightness profile as 
S oc , then typical globular clusters have slopes of 
77 ~ 4 (jMcLaughlin fc v an dcr Marcl 2005). We measure 
a shallow slope for Pal 13 of?7 = 2.8±0.3. This is con- 
sistent with our inferred King tidal radius which is two 
times larger than the calculated tidal radius of Pal 13 
based on its stellar mass and distance from the Milky 
Way (see §[3|). Our inferred slope is steeper than the C02 
estimate of ?7 = 1.8±0.2. These two profiles are compared 
in Figure |3l A key difference between the two King pro- 
files is the CMD region in which they were determined. 
While C02 determined the profile using stars inside a 
box in CMD space, we have used a window around the 
best-fitting isochrone. At faint magnitudes, there is sig- 
nificant contamination from unresolved galaxies, and our 
CMD window should reduce this contamination. 

In the right panel of Figure [21 we show the spatial 
density map of Pal 13. The outer isophotes in this figure 
are 3cr contours above the background, corresponding to 
r ~30 mag arcsec"^. The red dashed line/circle in both 



panels of Figure [31 are drawn at the King tidal radius rt — 
13.9' (98 pc). The photometric isophotes are regular out 
to the half-light radius of the cluster, and appear mildly 
irregular outside this radius. This is in contrast to other 
faint globular clusters such as Palomar 1 and Palomar 5 
which show evidence for tidal tails at much brighter 
surface brightness levels ([Niederste-Ostholt et al.l [20101 : 
lOdenkirchen et all 120031) . although similar to that ob- 
served for Palomar 14 (jSollima et a"L[|2011[ ). We inter- 
pret the slightly irregular isophotes and shallow density 
profile as mild evidence for tidal disruption and discuss 
further in §[S1 

Finally, we estimated the total luminosity of Pal 13. 
The published valu es for Pal 13 range between 1.2 x 10 ^ 
and 3.5 x 10^ Lq (ISieeel et al.l [200l[: ICote et al.l[200l) . 
making traditional methods of adding stars' individual 
fluxes too sensitive to the inclusion (or exclusion) of po- 
tential members (outliers). To alleviate these shot noise 
issues, the method used here relies solely on the total 
number of stars belonging t o the satellite and not on 
their individual magnitudes (jMunoz et al.ll2010l) . To es- 
timate an absolute magnitude, we modelled the satel- 
lite's population with the best-fitting theoretical lumi- 
nosity functio n, in this case a 12 Gyr population with 
[Fe/H]= -1.6 ([Girardi et al.|[2002D . We then integrated 
the theoretical luminosity function to obtain the total 
flux down to a given magnitude limit. The last step was 
to scale this flux using the total number of Pal 13 stars in 
our catalog down to the same magnitude limit. To esti- 
mate the uncertainty in the absolute magnitude we car- 
ried out a bootstrap analysis generating 10,000 realiza- 
tions of Pal 13 from the ph otometric data (for more de- 
tails see lMunoz et al.|[2010D . We assu med a Chabrier ini- 
tial m ass function for the estimation (iBruzual fc CharlotI 
120031 ). We obtained My = -2.8 ± 0.4, or equivalently 
Lv = 1.1^0 3 ^ ^0. The central surface brightness is 
^ay — 23.9lQ ymag arcs~^. While our results are one 
sigma lower than published values, it is a more robust 
estimate of the true luminosity of Pal 13 as it better ac- 
counts for shot noise in this quantity. We compare the 
size and luminosity of Pal 13 to other Milky Way stellar 
systems in Figure [1] 

3. MEMBERSHIP CRITERIA 

The systemic velocity of Pal 13 lies within the velocity 
distribution of foreground Milky Way stars. We therefore 
selected spectroscopic members of Pal 13 by combining 
color-magnitude and velocity criteria with proper motion 
information from the literature. 

We flrst applied a color-magnitude selection window 
shown in Figure [H CMD selection was perfor med using 
a 12 Gyr, [Fe/H] = -1.6 lGirardi et al.l (|2002[ ) isochrone 
in the CFHT photometric system, shown in Figure [21 to 
flt Pal 13. The selection window was defined as 0.1 mag 
around this isochrone added in quadrature to the 2-a 
photometric errors. The CMD selection window does not 
inclu de blue strag gler stars which are abundant in Pal 13 
(Clar k et al.l 12004^ . We selected blue stragglers sepa- 
rately by defining a generous window bluer and brighter 
than the main sequence turnoff, but dimmer than the 
horizontal branch (Figure [3]). Stars in the blue straggler 
window are removed from all analysis below unless spec- 
ified. We study the blue straggler population separately 



5 




10 10.0 OA 0.3 0.2 0.1 0.0 -0.1 -0.3 -0.3 -0.4 

radius (arcmin) Ara (°) 



Fig. 3. — (Left) Pal 13 surface brightness profile based on CFHT/Megacam imaging reaching a magnitude limit of r = 24. Here we 
compare our best fitting King profile (solid) to that from C02 (dotted) and parameterize the outer slope of the surface brightness profile 
(dashed). The dot-dashed line represents the background level. (Right) The surface density profile for Pall3. The outermost contours 
are our 3cr confidence limits, corresponding to 30 mag arcsec"^. The dashed red line in both panels is plotted at the King tidal radius of 
rt = 13.9' = 98 pc. 
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Fig. 4. — (Left) CMD of Pal 13 based on CFHT imaging for all spectroscopic stellar targets. The primary color selection window is 
plotted in light red with the best-fitting isochrone centered as a black dotted line. The blue straggler selection window is plotted in light 
blue. Red circles represent likely velocity members within the primary selection window, red squares in this region are further confirmed 
as Pal 13 members via proper motions. Small black dots represent definite non-members. We also identify two stars from the C02 study as 
orange stars. ORS-91 lies near the horizontal branch and ORS-38 lies near the RGB. (Right) Spatial diagram of Pal 13 for all spectroscopic 
stellar targets. Symbols are the same as the CMD plot on left. The solid circle is the cluster's half-light radius (ri^2 = 0.9'). The dashed 
red circle represents the King tidal radius (rt = 13.9'). 



in S OI 

After the CMD selection, we applied a basic velocity 
cut to the sample. We chose a conservative velocity win- 
dow of 10-a around the systemic velocity as measured by 
C02. We use the velocity dispersion measured by C02 of 
(7 = 2.2 km s~^, thus our velocity window corresponds to 
±22 km around the systemic velocity of 26.2 km s~^. 
Combined with stars passing the CMD cut, this leaves 



93 member candidates of Pal 13. 

Milky Way foreground stars are expected to contami- 
nate our CMD and velocity selected sample. Using th e 
Besangon model of the Milky Way (|Robin et al.ll2003| ). 
we applied the same selection criteria and predict 5 to 
7 foreground stars in our sample of 93 Pal 13 members, 
depending on details of normalizing the model to data. 
A maximum of 6 of these possible 7 foreground stars are 
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predicted in the thin and thick disk, and 1 in the halo. 
A method to identify foreground stars is via metahicity, 
since the average [Fe/H] of Pall3 is significantly lower 
than expected for the majority of foreground Milky Way 
disk stars f§ l4.4|) . As explained in § 14.41 we can measure 
metallicities for only the 20 brightest stars in our sam- 
ple. As seen in Figure [71 four stars are far more metal- 
rich compared to the main sample. Three stars are at 
very large radius from Pal 13 where the contamination is 
greatest. We reject these four stars as members, leaving 
89 candidate PallB stars in our 'full' sample. 

Another way to identify foreground stars is via proper 
motions. Proper motions are ava ilable for 69 out o f 
89 candidate members (78%) from iSiegel et al.l (|2001[ ). 
Three stars in this sample have large proper motions that 
are inconsistent with being cluster members (assigned 
0% membership probabilities by Siegel et al.), and we 
eliminate them from our sample. This leaves 66 can- 
didate member stars of Pal 13. We defined this sample 
of 66 stars as our 'cleaned' Pal 13 sample, and proceed 
with analysis on these stars. We cannot be certain the 
cleaned sample is free of Milky Way contamination, but 
estimate from the Besangon models that there is likely 
less than one foreground star given our sample criteria. 
The results of this selection are shown in Figure [SJ 

Tidal heating or evaporation processes will unbind 
stars from the gravitational potential of Pal 13. These 
stars can be a source of contamination if they remain 
in the spatial vicinity of the clus ter and would not 
be re jected by our methods above |Klimentowski et al.l 
|2007[) . While the three stars at largest radius passing 
our CMD and velocity cuts are prime candidates for un- 
bound Pal 13 stars, these were rejected based on their 
high metallicity and are more likely to be Milky Way 
foreground stars. 

To estimate at what projected radius stars are likely to 
be unbound from Pal 13, we calculate the Jacobi tidal ra- 
dius assuming a Galactocentric distance of 25.3 kpc and 
a Milky Way mass inside this radius of 3 x 10^^ Mq. 
The Jacobi tidal radius of Pal 13 is 6.7' (47 pc), larger 
than the half-light radius of the cluster rcff = 1.3', but 
smaller than the inferred King tidal radius of/'t = 13.9', 
as noted pre v iously by C02, ISiegel et"aLl ()2001[ ) and 
iKiipper et al.l ()2011h . Although the Jacobi radius pro- 
vides only a very approximate estimate of the tidal ra- 
dius (jBinnev fc Tremainell2008l ). this suggests that stars 
in the outskirts of Pal 13 may not be bound to the clus- 
ter. This calculation assumes a circular orbit. Stars in 
our 'clean' kinematic sample are located well inside the 
Jacobi tidal radius, while our full sample contains two 
stars beyond this distance. 

3.1. Repeated Velocity Measurements 

We obtained numerous repeated velocity measure- 
ments of stars in Pal 13 in order to investigate the fre- 
quency of unresolved binary stars and their influence on 
the velocity dispersion. We measured two or more veloc- 
ities for 62 out of 306 total (36 out of 89 member) stars. 
We further fold in historical data from C02, taken ten 
years prior to our DEIMOS observations, providing mul- 
tiple measurements for 72 out of 306 total (42 out of 89 
member) stars. Binary fractions inferred from our data 
are necessarily lower limits on the true binary fraction 
and are most sensitive towards short (few day) periods. 



but also include one year to ten year periods depending 
on the available data. We also note that the repeated 
stars are heavily biased towards bright magnitudes. 

Using the error-normalized velocity difference plotted 
in Figureini we define a star to be variable if individual in- 
dependent measurements are more than 3cr discrepant in 
this quantity. We acknowledge that the number of 'iden- 
tified' binaries depends sensitively on this choice. We 
find that 5 out of 89 in our full member sample are ve- 
locity variables (6%). These 5 stars are also contained in 
our cleaned sample. Given the large velocity differences 
(a few to tens of km s~^), we assume these systems are 
unresolved binary stars, as single variable stars are un- 
likely to be in this region of color-magnitude space. As 
expected, the fraction of binary stars identified in our 
incomplete s ample is far less than 30 ± 4% measured by 
IClark et al.l (j2004[) . but does provide spectroscopic con- 
firmation that unresolved binary stars exist in Pal 13. 

For all stars with multiple measurements, we use the 
error-weighted mean velocity in the calculations below. 
Because we are measuring the five binary stars at a ran- 
dom points in their orbits, it is not guaranteed that this 
will be the true systemic velocity of the star. We will 
therefore calculate, in § 14.21 the velocity dispersion of 
Pal 13 with and without the 5 identified binary systems. 

4. RESULTS 
4.1. Comparison to Cote et al. (2002) results 

The velocity dispersion of Pal 13 measured by C02 of 
2.2±0.4 km s~^ is significantly larger than that predicted 
based on the stellar mass alone of 0.3 — 0.4km s~^. 
If this velocity dispersion accurately traces the motion 
of stars in a gravitational potential, it implies over 
95% of mass in Pal 13 is not luminous, i.e. that Pal 13 
has a significant dark matter component. We mea- 
sured Keck/DEIMOS velocities for all 21 stars in the 
Keck/HIRES C02 sample and, in most cases, have mul- 
tiple measurements separated by up to one year. 

Throughout this paper, we determine dispersions us- 
ing a Monte Carlo Markov C hain (MCMC) method 
(|Lane et al.ll2010l : lGregorvll2005f ) assuming uniform pri- 
ors on the velocity and velocity dispersion. This method 
produces the same res ults as the likel i hood maximization 
technique described in lWalker et al.l (|2006[ ) when the er- 
ror distributions are Gaussian. In the regime where the 
velocity dispersion approaches the observational errors 
on this quantity, the MCMC method allows us to calcu- 
late confidence intervals directly from posterior distribu- 
tions without the assumption of Gaussianity. We explore 
in greater detail the choice of priors and sampling meth- 
ods in this regime in Geha et al. (in prep). We confirmed 
that our MCMC algorithm reproduces the same disper- 
sion and errors reported by C02 for their HIRES dataset. 

Using a single epoch of velocities and the 21 star sam- 
ple defined in C02, we measured a Keck/DEIMOS veloc- 
ity dispersion of 2.5±0.7km s"^, in good agreement with 
the Keck/HIRES measurement. We next combined our 
47 independent measurements of these 21 stars (six stars 
have only a single measured velocity, while 3 stars have 
four or more epochs), and redetermine the velocity dis- 
persion. Combining our repeat measurements, we find a 
lower velocity dispersion of 1.6 ± 0.7 km s""'^, suggesting 
that unresolved binary stars are indeed inflating the sin- 
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Fig. 5. — (Left) Velocity histogram of Pal 13 stars passing our primary color-magnitude selection criterion. We use a bin size of 2.5 km s"^ 
representing the typical velocity error. The normalized Milky Way model is plotted as a solid blue line and estimates the foreground 
contamination in each bin. CMD selected non-members are plotted as a solid black line and excludes blue stragglers. (Right) Velocity 
of CMD-selected stars as a function of cluster radius. The dashed line represents the half-light radius. Red circles represent CMD and 
velocity selected Pal 13 stars, red squares are stars which are further confirmed as Pal 13 members via proper motions. Small black dots 
show non-members. We also identify binaries as blue diamonds. 



gle epoch dispersion, although not to a level consistent 
with the stellar mass alone. 

We next examined the C02 stars for evidence of vari- 
ability. There are two stars in this sample which we 
suspect, ORS-38 and ORS-91, identified as orange sym- 
bols in Figure m ORS-38 is near the RGB and varies in 
velocity by 17 km s~i between the DEIMOS and HIRES 
measurements. This is one of the five variables identi- 
fied in §0 While it does not inflate the HIRES ve- 
locity dispersion, ORS-38 is the primary star inflating 
our calculations. ORS-91 is on the horizo ntal branch 
and li es in the instability band identified in ' Ivezic et alj 
(|2005[) . Ahhough the DEIMOS and HIRES velocities 
agree, and including/removing this stars does not affect 
our results, we chose to exclude this and other horizon- 
tal branch stars from our sample. Removing ORS-38 
and ORS-91 does not affect the HIRES-determined dis- 
persion, but the DEIMOS-based velocity dispersion falls 
to 0.7lQg km s~^. This dispersion is consistent with 
a normal stellar mass-to-light ratio, although our mea- 
surement uncertainties are too large to rule out a higher 
than expected mass-to-light ratio. We plot the error- 
normalized differences in velocity between the HIRES 
and DEIMOS samples in Figure [6l noting the position 
of these two rejected stars. 

4.2. The Velocity Dispersion of Palomar 13 

We determine the velocity dispersion of Pal 13 first us- 
ing the cleaned sample of 66 stars defined above. This 
sample has minimal Milky Way foreground contamina- 
tion and contains 5 known velocity variable stars. Using 
only a single epoch of DEIMOS observations, we calcu- 
late a dispersion of 2.3 ± 0.5km s^^. Using a weighted 
average of all our repeat measurements, the velocity dis- 
persion decreases to 1.6 ± 0.5 km s~^. By removing the 
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Fig. 6. — Data comparison to Cote et al. HIRES velocities. We 
plot the histogram for the velocity difference, normalized by obser- 
vational error: (■yDEIMOS-'!'HIREs)(o"DEiM0S+'^HIREg)~^^^- "^^^ 
two Cote et al. objects in question are identified in their respective 
bins (ORS-91 and ORS-38). 



5 known variable stars, the velocity dispersion drops 
close to our measurable limits: a = 0.6^05 km s~^. 
Finally, adding in quadrature measurements from C02 
where available, the velocity dispersion of Pal 13 is cr = 
These results demonstrate that unre- 



0.4l° jkm s- 



solved binaries or variables, even multiply measured and 
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averaged, can still contribute significantly to the disper- 
sion given our sample size. 

We take as the final velocity dispersion, the value de- 
termine without binaries combining both DEIMOS and 
HIRES measurements (61 stars, a ~ 0.4lQ'4km s~^). 
The average cluster distance of this sample is 1.5', 
slightly larger than the Pal 13 half-light radius of ri/2 = 
1.3', as seen in right panel of Figure [51 We do not have 
a sufhcient number of stars to produce a binned velocity 
dispersion as a function of radius, however, splitting the 
sample at the radius enclosing half the measured stars, 
we find that both the inner and outer bin have similar 
velocity dispersions, but the error on each value doubles. 

4.3. A Stellar Mass-To-Light Ratio 

We calculate the mass and mass-to-light ratio of 
Pal 13 using the velocity dispersion determined via our 
61 member clean sample without binary stars, a = 
0.4lQ4km s~^. We determine the dynamical mass en- 
close d within the half-light radius applying the formula 
from IWolf et all (poToh: Mi/2(ri/2) ^ AG-\al^)n,2, 
where M1/2 is the mass contained within the 3D pro- 
jected half-light radius, a\os is the line-of-sight veloc- 
ity dispersion. This mass estimator is based on the 
Jeans Equation and is less sensitive to uncertainties 
in the velocity anisotropy than other simple mass es- 
timates. We calculated an enclosed mass of Mi/ 2 < 
1.3li'3 ^ ^0- Using the luminosity determined in 
ii l2.2.21 Lv = l-lio;3 ^ 10^ Lq, we calculated a mass-to- 
light ratio of Ty = 2.412;° M©/ L© within the half-light 
radius. 

We alternatively calculate the mass-to-light ratio of 
Pal 13 assuming velocity isotropy and that the clus- 
ter is well fit by a King profile. Following C02, the 
mass-to-light ratio is given b y: Ty = aQ/{0.003rcfJ,o,v) 
(jRichstone fc Tremaind[T986[ ). where Tc is the King core 
radius in parsecs, /io,v the central surface brightness 
( L0 pc^^) and (To the central velocity dispersion. Using 
the best-fitting King parameters from ij |2.2.2[ we deter- 
mine Ty = 1.9 [need to calc errors]. 

Within the limits of our velocity precision, the mass- 
to-light ratio of Pal 13 is consistent with its stellar mass. 
Our estimates are consistent with the average dynamical 
mass-to- ligh t ratio of old Milky Way globular clusters of 
Ty = 1.45 (IMcLaughlinll200fli) . Since Pal 13 is slightly 
more metal-poor than the average Milky Way globular 
cluster in this sample, we computed the theoretical mass- 
to-light ratio, based on a Chabrier in itial mass function, 
an ag e of 12 Gyr and [Fe/H] = -1.6 (jBruzual fc CharlotI 
I2OO3I) to be Ty = 2.0. We therefore conclude that the 
dynamics of Pal 13 are consistent with its stellar mass 
alone 

4.4. Spectroscopic Metallicity 

We measured spectroscopic metallicities for individual 
stars in Pal 13 using the spectral synthesis method intro- 
duce d bv lKirbv et al . (2008a) and refined by Kirb v et al.l 
(|201C?) . Metallicities based on this method are reliable for 
stars with log 5 < 3.6 (at the distance of Pal 13 roughly 
r < 20). In Figure [71 we plot metallicities of the 20 
stars passing this criterion in our full sample. As dis- 
cussed in §[3l four stars have clearly deviant metallicities 
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Fig. 7. — Spectroscopic metallicity as a function of cluster dis- 
tance. We reject as Pal 13 members the four stars with metal- 
rich measurements. The average metallicity of Pal 13 is [Fc/H] 
= —1.6 ± 0.1, with an internal metallicity spread 0.1 ± 0.1 dex. 

from the main sample and we have excluded these stars 
from Pal 13 membership. The weighted average metallic- 
ity for the 16 member stars is ([Fe/H]) = —1.6 ±0.1 dex. 
This is comparable to previous spectroscopic estimates of 
the metallicity based on medium resolution s pectroscopy, 
[Fe/H] = -1.67 ±0.15 (IZinn fc piaz| [T98l and [Fe/H] 
= -1.9 ± 0.1, (|Friel et al.l 119821 ). as well as [Fe/H] = 
— 1.98 ± 0.31 dex b ased on a single Pal 13 star observed 
at high resolution (jCote et al.ll2002D . 

Using the same sample of 16 stars above, we esti- 
mated an internal metallicity dispersion of 0.1 ± 0.1 dex. 
While the one-sigma errors are consistent with no inter- 
nal metallicity dispersion, as expected for a single stellar 
population globular cluster, the measurement does al- 
low a modest metallicity dispersion. This dispersion is 
smaller than obs erved in some lumin ous globular clusters 
(e.g., NGC 2419 ICohen eFall I20T0I ) . and is significantly 
smaller than observed for stellar systems containing dark 
matter, e.g. dwarf galaxies, sh ow clear metallicit y disper- 
sions on the order of 0.5 dex (jKirbv et al.ll2010D . Dwarf 
galaxies also exhibit a tight correlation b etween luminos- 
ity and metallicity ([Kirbv et al.ll2008bD . Given Pall3's 
low luminosity (My = —2.8), this relationship predicts a 
metallicity of [Fe/H]^ —3, significantly more metal-poor 
than observed. Both the average metallicity of Pal 13 
and the low internal metallicity dispersion are further 
evidence that Pal 13 is a globular cluster. 

4.5. Blue Straggler Stars 

We show above that the dynamics of Pal 13 are in- 
fluenced by the presence of unresolved binary stars. 
This is unlike other Milky Way globular clusters whose 
single-epoch velocity dispersions arc consistent with their 
stellar mass ([B aumgardt et al. 2009; Lane et al. 201(J 
IHankev fc GaS l2011f ). We next investigate the blue 
straggler stars (BSSs) population of Pal 13 as a means 
to determine whether the binary fraction in Pal 13 is 
indeed higher than typical globular clusters. BSSs are 
positioned brighter and blueward of the primary main- 
sequence turnoff. In globular clusters, BSSs mimic a 
younger stellar population, but are interpreted as sys- 
tems which formed via stellar collisions (Hills & Day 



1976 ) or mass transfer in close binary systems (McGreal 
1964f ). We have not included BSSs in our kinematic 



analysis, however, our spectroscopic observations of these 
stars provide some clue to their formation, and their fre- 
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quency provides further evidence for an enhanced binary 
fraction in Pal 13. 

In our spectroscopic sample, 12 stars pass the BSS cri- 
terion in color-magnitude (Figure |4l). Two of these stars 
do not pass our velocity criterion for membership, and 
their proper motions measured by Siegel et aL (2001,) 
confirm they are not Pal 13 members. The remaining ten 
stars have proper motions consistent with Pal 13 mem- 
bership. We obtained repeated velocity measurements 
for only one member star; the repeat measurements are 
consistent with a constant velocity. However, the veloc- 
ity distribution of the 10 stars identified as BSS members 
of Pal 13 is broader than the overall velocity distribution 
of Pal 13. This is reflected in the velocity dispersion of 
the BSSs, 4.1 ± 1.4 km s~^ substantially larger than the 
main Pal 13 sample. This hints that at least some of 
these BSSs may be close binary systems. 

A common way to express the BSS abundance is by 
measuring the ratio, of BSSs to RGB stars, F'^^l. Using 
our CFHT photometric data, we estimated the fre quency 
of BSSs following the method of Leig h et al.l 
(|2007f ) by counting the number of BSSs and RGB stars 
within well defined regions of the CMD. These regions 
are similar to Leigh & Sills, with the main difference be- 
ing our BSS box is located closer to the main sequence 
turnoff. The reason for this is that Pal 13 is a low density 
cluster and overcrowding near the center is not an issue 
and, therefore, the scatter in the main sequence is lower. 

We measured F'^'^Jl in both a region within the half- 
light radius (1.3') and three times this value obtaining 
F^^l = 1.0 ± 0.4 and 0.8 ± 0.3, respectively. Our resuhs 
are directly comparable to the values of determined 
by iLeigh et al.l (|2007[) for 57 Milky Way globular clus- 
ters (their sample did not include Pal 13). These authors 
found an average value of F'^^l — 0.14 with a standard 
deviation of 0.19. In comparison. Pal 13 presents a re- 
markably high fr equency of BSSs. G iven the binary frac- 
tion measured bv lClark et"all (|2004l ) of 30 ±4% based on 
stars redwards of the main sequence, the high frequency 
of BSS stars can be interpreted as due to mass transfer 
between primordial binaries. The high BSS frequency 
is consistent with Pal 13 having a higher binary fraction 
than typical Milky Way globulars. 

5. DISCUSSION AND CONCLUSIONS 

We demonstrate that the mass-to-light ratio of the 
Milky Way globular cluster Pal 13 is consistent with 
its stellar mass, based on deep CFHT/MegaCam imag- 
ing and Keck/DEIMOS spectroscopy. Our spectroscopic 
member sample is triple that of previous studies, consist- 
ing of 61 stars with multiply measured velocities on the 
timescale of days and years. We have minimized contam- 
ination from Milky Wa y foregro und stars via proper mo- 
tion estimates from . Siegel et al.l (|2001), and removed un- 
resolved binary stars when possible via multiple spectro- 
scopic measurements. When combined with HIRES ve- 
locities from C02, our final velocity dispersion for Pal 13 
is (T = 0.4lg3km s~^, corresponding to a mass-to-light 
ratio of Ty = 2.4t2'4M0/Lo. Thus, the dynamical mass 
of Pal 13 is consistent with its stellar mass and there is 
no need for dark matter. We t herefore confirm the sug- 
gestion by ,Blecha et all (|2004f) that binaries inflate the 



single epoch velocity dispersion. While Pal 13 appears to 
be a 'normal' globular cluster in this sense, two indepen- 
dent observed quantities in Pal 13 are atypical: (1) unlike 
other Milky Way globular clusters, the single epoch ve- 
locity dispersion is higher than predicted based on the 
stellar mass alone, and (2) the surface brightness slope 
is shallower than typical globular clusters. We examine 
each of these statements below. 

The single epoch velocity dispersion of Pal 13 is in- 
flated due to unresolved binary stars in both the HIRES 
data from C02 and our own DEIMOS dataset. We ar- 
rive at our final dispersion by removing these binary 
systems. Yet unresolved binaries do not affect the ve- 
locity dispersion of oth e r Milky Way glo bular clusters 
(|Baumgardt et al.l l2009t iLane et all I2010D , at least not 
to the same extent. Thus, Pal 13 must have a higher bi- 
nary fraction as compared to these globular clusters. Es- 
timates from the main sequence (Clark et al...2004) and 
blue straggler stars (§|13] Santana et al. in prep) are con- 
sistent with this statement. However, we are necessarily 
comparing Pal 13's binary fraction to luminous and more 
concentrated globular clusters where such estimates ex- 
ist. 

A stellar mass-to-light ratio for Pal 13 does not pro- 
vide an explanation for the second anomalous obser- 
vations in Pal 13, the observed shallow surface bright- 
ness (77 = —2.8 ± 0.3), nor the large King tidal radius 
{rt = 98 ± 11 pc) which is significantly larger than the 
calculated Jacobi tidal radius of 47 pc. The half-light ra- 
dius of Pall3 (ri/2 = 1.27±0.16' = 9.0±l.lpc) is larger 
than the average Milky Way cluster, as well as the me- 
dian size of globul ar clusters in external galaxies (e.g., 
[Jordan et alTl200l [Masters et al.|[2010f ). Pall3 may be 
a Galactic analog of the "di ffuse star cluster s" seen in 
some other galaxies (jBrodie & Larscn 2002; Peng et aD 
l200fih . 

The photometric isophotes of Pal 13 (Figure El show 
some irregularity near the Jacobi radius, but we do not 
observe the well defined S-shape curve seen in other glob- 
ular clusters which are actively being tidally stripped. 
While Pal 13 may not currently be undergoing tidal heat- 
ing (e.g. lack of S-curve in photometry, low velocity dis- 
persion), tidal debris in the vicinity of the cl uster may 
influence its appearance. iKiipper et al.l (|2011[ ) proposed 
that Pal 13 is at apogalacticon (farthest approach), and 
the photometric properties can be explained by tidal de- 
bris stripped from the cluster throughout its orbit and 
compressed. Conflrming whether the photometric fea- 
tures of Pal 13 are intrinsic or due to tidal interactions 
may be tested by future studies that consist of measuring 
velocities and metallicities of stars at larger radius then 
the present sample, particularly between the predicted 
and King tidal radius (6 — 13') of Pal 13. 

While the large inferred radii of Pal 13 may be due 
to tidal interactions, its present luminosity and size are 
currently more closely related to the dark matter dom- 
inated Milky Way ultra- faint galaxies (Figure [Ij . Is it 
possible that the ultra-faint galaxies are also affected 
by tidal interactions or binary stars, and that these 
objects may instead hav e mass-to- light more simila r 
to globular clusters (e.g., iMcConnachie fc Cotel l2010l) ? 
This exercis e has been done f or the ultra- faint g alax- 
ies Segu e 1 ([G eha et al.' '2001 iSimon et al.l [20Tl[) and 
Bootes I (|Koposov et al.. ,201 ij) . whose velocities disper- 
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sions arc between 2-4 km s^^. In both cases, the binary- 
corrected velocity dispersions are still consistent with 
a significant dark matter content. Furthermore, both 
galaxies have a more significant spread in metallicities as 
compared to Pal 13. Thus, both the corrected velocity 
dispersions and metallicities suggest that the ultra-faint 
galaxies are distinct from low luminosity globular clus- 
ters such as Pal 13. 

The dynamical mass of Pal 13 inferred from the 
binary-corrected velocity dispersion is consistent with 
its stellar mass alone. However, the large influence of 
unresolved binaries on the velocity dispersion, as well 
as high inferred blue straggler population makes Pal 13 
a unique system. In addition Pal 13 has an unusual 
density structure suggesting that tidal interaction may 
be important in this low-mass halo object. Studying 
other low luminosity globular clusters in the Milky Way 



halo and beyond will answer whether these properties 
are miique to Pal 13 or are typical properties of low 
luminosity globular clusters. 
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TABLE 1 

Keck/DEIMOS Multi-Slitmask Observing Parameters 



Mask 


a (J2000) 
(h : m : s) 


(5 (J2000) 
(°:':") 


Date Observed 


PA 
(dog) 


^cxp 


# of slits 


% useful 
spectra 


PallS.l 


23:07:03.4 


+12:50:55 


Sept 4, 2008 


38 


9000 


74 


86% 


Pall3_2 


23:06:50.8 


+ 12:49:09 


Sept 4, 2008 


-13 


7200 


81 


95% 


Pall3.3 


23:06:47.3 


+ 12:44:14 


Sept 5, 2008 


112 


7200 


72 


94% 


Pall3.4 


23:06:41.7 


+12:43:45 


Sept 5, 2008 


180 


8100 


83 


95% 


Pall3_5 


23:06:08.4 


+ 12:49:16 


Sept 4, 2008 


17 


2700 


52 


96% 


Pall3_6 


23:07:22.7 


+12:49:25 


Sept 4, 2008 


-166 


2460 


56 


98% 


Pall3.7 


23:06:31.2 


+12:37:44 


Sept 5, 2008 


87 


2400 


51 


92% 


Pall3.8 


23:06:42.3 


+ 12:45:54 


Sept 5, 2008 


74 


1680 


42 


83% 


paI13 


23:06:51.0 


+12:45:14 


Oct 13, 2009 





1800 


23 


96% 


Note. — 


The right ascension, declination, date of observation. 


position 


angle and total expo- 



sure time for each Keck/DEIMOS slitmask. The final two columns refer to the total number of 
slitlets on each mask and the percentage of those slitlets for which a redsliift was measured. 



TABLE 2 

PALOMAR 13 PROPERTIES 



Property 



Cote ct al. (2002) 



This work 



Core Radius (King Model), 
Tidal Radius (King Model), rt 
Half-Light Radius (King Model), 
Luminosity, Ly 
Absolute Magnitude, My 
Metallicity, [Fe/H] 

Right Ascension (Center, King Model), 
Declination (Center, King Model), dec 
Distance Modulus, (m — M)q 
Distance, d 



0.23 + 0.03 arcmin 
26 + 6 arcmin 
N/A 

(2.8 + 0.4) X 10=^ Lq 
—3.8 mag 
-1.9 + 0.2 dex 
346.6867 deg 
12.7717 deg 
16.93 + 0.10 mag 
,+1.2 
-1.1 



24.3^ 



kpc 



0.42 + 0.06 arcmin 
13.9 + 1.50 arcmin 
1.27 + 0.16 arcmin 

(l.ltO j) X 103 Lq 

-2.8 + 0.4 mag 

-1.6 + 0.1 dex 

346.68519 + 0.00063 deg 

12.771539 + 0.00068 deg 

N/A 

N/A 



Note. — Wc list properties of Pal 13 compiled by C02 and as independently measured by this 
work for comparison. 



TABLE 3 

Keck/DEIMOS Velocity Measurements for Palomar 13 



i 


Name 


a (J2000) 
(h m s) 


(5 (J2000) 

(O , „) 


r— mag 
(mag) 


(s - 1') 

(mag) 


^avg 

( km s" 




( km s" 




( km s^^) 


Full 


Clean 





LRIS.61 


23:06:43.15 


+ 12: 


12:50.30 


19.85 


0.44 


29.85 + 


2.85 


29.85 + 


2.85 




N 


N 


1 


LRIS.21 


23:06:42.00 


+ 12: 


12:26.50 


17.68 


0.75 


20.27 + 


0.26 


27.13 + 


2.25 


20.18 + 0.26 


Y 


Y 


2 


LRIS.47 


23:06:45.59 


+ 12: 


12:39.50 


19.36 


0.68 


23.37 + 


0.69 


35.15 + 


2.23 


22.14 + 0.72 


Y 


Y 


3 


LRIS.22 


23:06:48.27 


+ 12: 


12:46.60 


18.30 


0.61 


19.98 + 


0.46 


25.27 + 


2.21 


19.74 + 0.47 


Y 


Y 




















19.61 + 


2.23 


19.93 + 0.54 
























29.36 + 


2.22 


19.41 + 0.93 
























25.72 + 


2.22 


18.87 + 1.84 
























21.50 + 


2.27 








4 


LRIS.64 


23:06:41.61 


+ 12: 


12:09.50 


19.91 


0.54 


30.42 + 


2.61 


30.42 + 


2.61 




Y 


Y 




















32.39 + 


2.85 


























27.43 + 


3.13 








5 


LRIS.41 


23:06:44.76 


+ 12: 


12:18.80 


19.75 


-0.12 


25.94 + 


3.41 


25.94 + 


3.41 




N 


N 




















25.94 + 


3.41 








6 


LRIS.15 


23:06:48.52 


+ 12: 


12:19.20 


17.32 


0.72 


29.22 + 


0.27 


24.68 + 


2.20 


29.29 + 0.27 


Y 


Y 




















23.92 + 


2.21 


29.14 + 0.80 
























25.01 + 


2.20 


30.55 + 0.64 
























24.65 + 


2.22 


28.76 + 0.36 




























29.80 + 0.62 






7 


LRIS.39 


23:06:43.96 


+ 12: 


12:20.20 


19.33 


0.58 


25.60 + 


0.68 


30.15 + 


2.38 


25.19 + 0.71 


Y 


N 




















30.15 + 


2.38 


25.19 + 0.71 







Note. 
apparent 



Velocity measurements, magnitude, color and membership data for spectroscopically observed stars of Pal 13. Position, 



-band magnitude, {g ■ 



color, weighted average of DEIMOS and HIRES heliocentric velocities (-Ua 



DEIMOS he- 



liocentric radial velocity HIRES heliocentric radial velocity (fR) from C02, member and clean member status for each star as 

determined in §[2] and §[3] Each sub-row lists the subsequent DEIMOS and HIRES individual epoch velocity measurements. Note 
that the HIRES velocities have been shifted by +0.5 km s~^for zero-point alignment to this study. This table will be published in 
its entirety in the electronic edition of the Astrophysical Journal. 



